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In the normal metal/ferromagnetic insulator bilayer (such as Pt/Y 3 Fe 50 i 2 ) and the normal 
metal/ferromagnetic metal/oxide trilayer (such as Pt/Co/AlOx) where spin injection and ejection 
are achieved by the spin Hall effect in the normal metal, we propose a minimal model based on 
quantum tunneling of spins to explain the spin-transfer torque and spin pumping caused by the spin 
Hall effect. The ratio of their damping-like to field-like component depends on the tunneling wave 
function that is strongly influenced by generic material properties such as interface s — d coupling, 
insulating gap, and layer thickness, yet the spin relaxation plays a minor role. The quantified result 
renders our minimal model an inexpensive tool for searching for appropriate materials. 

PACS numbers: 75.76.-l-j, 75.47.-m, 85.75.-d 


Introduction. - Two reciprocal mechanisms, namely the 
spin-transfer torque (STT)[1, 2] and spin pumping[3, 4], 
play the major roles in spintronic devices. Besides in the 
originally proposed metallic heterostructures, they also 
manifest in heterostructures involving insulators, such as 
the normal metal/ferromagnetic insulator (NM/FMI) bi¬ 
layer realized by Pt/Y 3 Fe 50 i 2 (Pt/YIG)[5], where STT 
can be used to excite magnons[6, 7]. Unlike in metal¬ 
lic heterostructures, the usual way of spin injection by 
using a spin-polarized charge current in the current- 
perpendicular-to-plane (CPP) geometry is difficult in 
NM/FMI because the insulating FMI impedes the charge 
current. Instead, the spin injection in this system is done 
by using the spin Hall effect (SHE) in the NM, which in 
the current-in-plane (CIP) geometry injects a pure spin 
current into the FMI to cause STT. In the reciprocal 
process, SHE converts the pure spin current generated 
from spin pumping into electric signals. Despite pointing 
at promising applications in magnetic memory devices, 
the microscopic mechanism concerning how the pure spin 
current tunnels into an insulator to cause STT and spin 
pumping remains unclear, especially if spin relaxation 
plays a crucial role as in metallic heterostructures [1]. 

Another system of similar kind yet involves more com¬ 
plexity is the normal metal/ferromagnetic metal/oxide 
(NM/FMM/oxide) trilayer, such as Pt/Co/AlOa, and 
similar ones[10-14]. It is generally suspected that the 
observed spin torque contains both SHE-STT[13] and 
the spin-orbit torque (SOT) [15-17, 19] that stems from 
the charge current flowing through the parity-breaking 
FMM. Concerning SOT alone, a field-like component is 
expected from the current induced spin density [21], yet 
scattering-related mechanisms[17~20], as well as intrinsic 
Berry curvature [22] have been shown to cause a damping¬ 
like torque of similar strength. On the other hand, the 
mechanism of SHE-STT is less explored, thus the relative 


weighting between SOT and SHE-STT remains unclear. 
Moreover, it is technically important to single out the 
contribution from SHE-STT alone and investigate if this 
part can be enhanced by any means. 

In this Letter, we present a minimal model for the STT 
and spin pumping in NM/FMI and NM/FMM/oxide 
based on quantum tunneling of spins. We point out 
that the damping-like and field-like component (defined 
in Eq. (5)) depend on the wave function tunneled into the 
ferromagnet that is strongly influenced by generic mate¬ 
rial properties such as interface s — d coupling, insulating 
gap, and thickness of the ferromagnet, but not crucially 
on spin relaxation. Using the two methods that originally 
predict the STT in metallic heterostructures, namely the 
Landau-Lifshitz (LL) dynamics]!] and angular momen¬ 
tum conservation[2], the spin mixing conductance[27] for 
STT is calculated. At present, it is of particular inter¬ 
est to find ways to enhance the spin mixing conductance, 
which may realize magnetic memory devices with more 
efficient magnetization switching and lower power con¬ 
sumption, and our minimal model serves as an inexpen¬ 
sive tool to guide the search for suitable materials. The 
mechanism of spin pumping is further clarified based on 
the adiabatic assumption, yielding Onsager relation ex¬ 
plicitly satisfied within this minimal model. The consis¬ 
tency with various experiments will be explained thor¬ 
oughly. 

NM/FMI bilayer.- We first address the quantum tun¬ 
neling of spins in the NM/EMI bilayer. Consider the 
bilayer shown in Fig. 1 (a) that contains two regions: 
(1) A NM at — oo < X < 0 described semiclasically by 
Hn = /2m — tJ/, where = F\pi//2 is the spin volt¬ 

age of spin cr = {t, 4-} at position x caused by SHE [24- 
26]. Eor an up spin incident from the left, the wave func- 
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tion at a; < 0 is 


macrospin S. From the LL dynamics, the s — d coupling 
Hsd = Fcr • S renders the STT[1], whose response in the 
damping-like and field-like direction define the spin mix¬ 
ing conductance[27] and Gi, respectively, 


where k^a = \/‘2m {e + ^j,Q)/h, and e is the Fermi energy. 
We consider a charge current j°y, so electrons moving 
in X direction has f || z because of the SFIE relation 
fcocrX II cr X y, and gives a positive spin voltage /xq || 
z at the interface such that fcot > koi- (2) A FMI at 
a; > 0 described by Hpi = I2m -I- Vq -b FS • cr, where 
Vb > e is the potential step. The S • cr term describes the 
s—d hybridization of the conduction electron spin and the 
magnetization S = ^ (sin 0 cos , sin 0 sin </?, cos 0). The 
evanescent wave function in the FMI is 
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after substituting Eq. (17) and |Ap = A^f|MoI/o^- 
Alternatively, one can extract STT from the spin 
current [2]. The spin current in the NM at position x 
is 

jx = [V’wO- (da^lpN) - {dxi^N) O'V'iv] ■ (6) 


where q± = y^2m (Vb ± r*S' — e)/^ > 0. We choose F < 0 
such that S tends to align with cr, and q- > q+. Since 
the effective mass m merely rescales various energy pa¬ 
rameters (see caption of Fig. 1), it is set to be the same 
in the formalism below for simplicity. 

Spin-transfer torque.- The amplitudes A E are 
solved by matching the wave function and its derivative 
at X = 0 , leaving only one free variable \A\^ that is at¬ 
tributed to the Fermi surface-averaged spin density at the 
interface |Ap = Nplfiol/a^, where Np is the density of 
states per at the Fermi surface, and a is the lattice con¬ 
stant. The spin of conduction electrons inside the FMI 
can be calculated from the wave function in Eq. (2). It 
is customary[ 1 ] to express them in the frame {x 2 ,y 2 ,Z 2 ) 
defined in Eig. 1 (a), where 2,2 || S, y 2 = Ao x S/sin0, 
and X 2 = — S X ^Ao x /sin0. The spinors in Eq. (2) 

are simply (1 0)^ and (0 1)^ in this frame. The STT can 
be calculated by the following procedure [1]. Introducing 
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the conduction electron spin tunneled into the FMI is 
(cr) = {il)pi\cr\ifpi )whose components are 


= -4^sin0e-(«++«-)“ (Re,Im) (n|+n^_) , 

I 4|2 

(a^=) = 2fA^cos0e-('?++«-)^ 

\leV 

X [(l»^t-rl^-+ (l’^i-rl^ + l»^l-l^)cos0] . (4) 


The total spin per cross section channel is denoted by 
(^) = dx((T). The magnetization within the 

range of (ct), about 27r/ (< 7 + + q-) ~ a, is treated as a 


Conservation of angular momentum requires that the 
spin current at the interface to be equal to the STT, 
{jo — joo) = of jo = t[ 2, 27, 28], which is indeed satis¬ 
fied in this approach. Consequently, one can use Eq. ( 6 ) 
to extract Gr,i [37] 

r5*A^F ^ ^ 2jo cosy; 2jgsin ip 

h I/iQ I sin 20 I/iQ I sin 20 

TSNp ^ _ josinip jpcosy? 
h * l/xoj sin0 l/xol sin0 

which gives the same G^y as in Eqs. (5). In short, ei¬ 
ther extracting from the LL dynamics [I] in Eqs. (5), or 
requiring angular momentum conservation[ 2 ] and using 
Eq. ( 6 ), one obtains the same STT. 

Gr and Gi are independent from the azimuthal angle 
< 7 ?, but have small 0 dependence through ye which can 
be considered as higher order contributions (At most few 
percent. Their ratio does not depend on 0). The experi¬ 
mentally measurable quantity Gr,i x e^ jh at 0 = O.Stt is 
plotted in Fig. I as a function of the s — d coupling and 
insulating gap relative to the Fermi energy. Their ratio 

G^ ^ /” {ay^)dx ^ kpi {q+-q-) 

G. /“(a-)dx k^,^ + qpq. ’ 

strongly depends on the wave function tunneled into the 
FMI and is therefore highly influenced by the insulat¬ 
ing gap and the s — d coupling. In the supplementary 
material [37], it is shown that the (anisotropic) spin relax¬ 
ation, simulated phenomenologically by an exponentially 
decay factor [I], is not crucial to determine the damping¬ 
like to field-like ratio in real materials. This is very dif¬ 
ferent from the metallic heterostructures in which spin 
relaxation is the origin of STT and gives predominately 
a damping-like torque[I, 2]. 
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For a reasonable value oi s — d coupling, a very small 
fiQ is already sufficient to create a damping-like torque 
that can overcome Gilbert damping and cause ferromag¬ 
netic resonance (FMR) [5], so the observed precession fre¬ 
quency (or the gyromagnetic ratio) is almost unaffected 
by the small field-like torque which can nevertheless be 
nonzero. A previous study shows that at the Ag/YIG 
interface, the enhanced local magnetic exchange field at 
the interface enhances Gr[29], which can be interpreted 
as enhanced s — d coupling in our calculation and is fairly 
consistent with the {g_,fcoi} > < 7 - 1 - situation. In the 
Stoner limit {F, |/.to|} {e, Vb — e}, the field-like com¬ 

ponent dominates, in good agreement with Ref.29. 



1 . 1.5 2 . 2.5 3 . 1 . 1.5 2 . 2.5 3 . 


(Vo-e)/€ (V|,-e)/€ 

FIG. 1: (color online) (a) Definition of the coordinates at 
= 0 and energy scales near the NM/FMI interface. Blue 
line indicates the potential profile seen by the condnction elec¬ 
trons. (b) The ratio of (c) damping-like to (d) field-like com¬ 
ponent of the spin mixing conductance at 6 = O.Stt plotted as 
a function of the insnlating gap {Vo — e) /e and interface s — d 
coupling —VS/e relative to the Fermi energy. The color scale 
is in units of ^ jha? which is abont ~ 

depending on the lattice constant a. The white regions in 
(b)'~(d) are where one potential barrier is lower than the 
Fermi energy To — e + FS' < e, hence irrelevant to the problem. 
If the effective mass in NM and in FMI are different, the axis 
labels of the plots are replaced by —VS/t —^ —mpTS/mN^ 
and {Vo — e) /e —>■ mp (Vb — e) jrapie. 


Spin pumping.- The spin pumping at the NM/FMI in¬ 
terface can be addressed by solving the Bloch equation 
of conduction electrons in the presence of a dynamical 
magnetization [4, 5] 

^ + dJ, = ^Sx{a)-Tsf (9) 

where is a spin relaxation term. Assuming adiabatic 


process, i.e., the magnetization dynamics \dS/dt\ is much 
slower than any frequency scale {e/h,Vo/h,T/h} in the 
problem, the wave functions in Eqs. (I) and (2), as well 
as (cr) in Eq. (17), remain valid. As detailed in the sup¬ 
plementary material[37], by considering the modification 
of (cr) after a small time lapse St due to the magneti¬ 
zation dynamics and taking St —^ 0, the right hand side 
of Eq. (9) vanishes. Upon integrating over x and using 
joo = 0, Eq. (9) gives 
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where Gr,i are the same as those obtained from Eqs. (5) 
and (7). Here /cq is interpreted as a proximity induced 
spin accumulation. Comparing Eqs. (5) and (10), evi¬ 
dently the Onsager relation[27], which dictates that the 
damping-like and field-like response for STT and for spin 
pumping must be equal, is satisfied. 

A previous calculation[5] gives generically Gi = 0 in 
spin pumping, which is however at odds with the first 
principle calculation of STT that shows both Gr and Gi 
are nonzero [29], as Onsager relation is not satisfied. Our 
calculation indicates that the material-dependent wave 
function is the key to resolve this discrepancy, and the 
spin pumping in general has both field-like and damping¬ 
like components. If the magnetization dynamics is a 
precession around a common axis, as that induced by 
EMR[5], the field-like component time-averages to zero 
while the damping-like component does not (see Eig. 1(b) 
of Ref. 37), so the observed dc spin current may be well 
explained by the damping-like component alone which 
however does not rule out the possibility of nonzero field¬ 
like component. Moreover, the field-like component also 
contributes to the ac measurements [30-32] which identify 
a much larger ac spin current than the dc one. Thus none 
of the known experiments contradicts the possibility of 
having a field-like component. 

NM/FMM/oxide trilayer. - The above analysis can also 
address the quantum effects in the NM/FMM/oxide 
trilayer, which may not be captured by diffusive 
approaches [33, 34]. Consider the trilayer shown in 
Fig. 2(a) that contains three regions: (1) A NM occupy¬ 
ing —oo < X < —b described by Eq. (1) with the replace¬ 
ment /Tq —>• (2) A FMM in —6 < a; < 0 described 

by HpM = p^/2m -|- FS ■ cr. It is unnecessary to include 
a Rashba term (p x Esoc) • cr because p || Esoq II x 
in this problem, where Esoq is the effective electric field 
due to the parity-breaking in x-direction[15-17, 19]. We 
consider a FMM much thinner than its spin relaxation 
length b Igf such that the wave function 
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FIG. 2: (color online) Same as Fig. (1), but for 

NM/FMM/oxide trilayer, (b) to (d) are plotted in terms of 
the s — d coupling relative to the Fermi energy —TS/e and the 
thickness of the FMM layer b/a. 


remains valid, where k± = \J 2m(e rS')/?i. (3) The ox¬ 
ide region at x > 0 described by Hq = j^m -b Vb0(x), 
where Vq > e is the potential step. The wave function is 

( 12 ) 


where q = ^2m{Vo — e)/h. The wave functions render 
zero spin current at the FMM/oxide interface jo = 0, but 
finite at the NM/FMM interface j-b ^ 0. After matching 
the wave functions, we introduce 
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the spin expectation values in the FMM are then 
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The Onsager relation is again satisfied. Alternatively, 
STT can be calculated by requiring angular momentum 
conservation^] aP {j-b — joa) = oPj-b = r and using 
Eqs. (6) and (7) with the replacement jo —!> j-b, which 
yields the same Gr,i- 

The numerical results for Gr,i at 0 = 0.37r are shown 
in Fig. 2, where certain oscillations with respect to the 
s — d coupling and the FMM layer thickness are evident, 
signaturing the effect of quantum interference, while Gr,i 
do not strongly depend on the insulating gap. This is 
in accordance with the measurement in Ta/CoFeB/MgO 
that shows hints for a varying GrjGi when changing 
FMM thickness [35], although one should keep in mind 
that we consider the NM thickness to be much larger 
than its spin relaxation length so the injected spin current 
saturates[12, 36], and hence our result does not depend 
on NM thickness unlike the experimental case. The abso¬ 
lute magnitude of ~ 10“^® ~ is close 

to that observed experimentally, thus the SHE-STT con¬ 
tribution relative to the SOT should not be overlooked. 

In summary, the generic material properties in het¬ 
erostructures involving insulators, such as s — d coupling, 
insulating gap, and thickness of the ferromagnet, are 
shown to strongly influence the tunneling wave function 
and subsequently the damping-like and field-like compo¬ 
nent of the STT and the spin pumping caused by SHE. 
Spin relaxation, on the other hand, plays a relatively 
minor role. The quantum effects of the wave function, 
such as quantum tunneling and quantum interference, 
are quantified in our minimal model that incorporates si¬ 
multaneously the LL dynamics, angular momentum con¬ 
servation, and Onsager relation, thus a convincing and 
inexpensive model to guide the search for appropriate 
materials. 

We thank P. W. Brouwer, P. Gambardella, P. Horsch, 
S. Maekawa, M. Mori, T. S. Nunner, J. Mendil, Y. 
Tserkovnyak, G. Vignale, and H.-H. Lin for stimulating 
discussions. 


Supplementary Material 


I. DETAIL OF SPIN-TRANSFER TORQUE 
CALCULATION 


The STT and spin pumping can be calculated 
from Eq. (5) with the replacement —>■ 

and using Eq. (14). The resulting Gr,i has 
small 6 dependence through Yg that can be attributed 
to higher harmonic terms allowed by symmetry [14]. 


A. NM/FMI bilayer 

By matching the wave function in Eq. (1) and (2) of 
the main text and their derivative at the interface x = 0, 
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the scattering coefficients in the NM/FMI bilayer read 
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(15) 


where na± and Fg are defined in Eq. (3) of the main text. 
Since the particle flux is zero /cot \A\^ — fcot \B\^— fcoj. \C\^ = 
0 , there is no charge current in this problem but only spin 
current, making it clear that the formalism describes the 
spin injection caused by the pure spin current in SHE. 

To extract the spin mixing conductance Gr,i from an¬ 
gular momentum conservation, we identify spin current 
at a; = 0 with the STT 
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where we have used S x x /xoj = |mo| S x x zj 

and S X /i,o = |/i.o| S x z since the spin voltage 
is defined to be polarized along /xq || z = ( 0 , 0 , 1 ), 
and the fact that the magnetization points at S = 
(cos0cos(/?,cos0sin(/?, sind). Solving for Gr,i in terms 
of each component of the spin current yields Eq. (7) of 
the main text. 

We now address the issue of spin relaxation, which in 
general can be anisotropic. Since we are not aware of any 
convincing calculation of spin relaxation for the evanes¬ 
cent wave function tunneled into the FMI, we model 
the spin relaxation phenomenologically by an exponen¬ 
tially decayed factor on the spin expectation value, al¬ 
though our argument below does not require an expo¬ 
nential form. Erom Eq. (4) in the main text, adding 
anisotropic spin relaxation length Ayj} into the spin 
expectation value yields 
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X (Re,Im) , (17) 


which gives Eq. (5) in the main text with corrections 
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(Im, Re) 
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(18) 


rendering the Gr and Gi that depend on Xy^ and , re¬ 
spectively, so the damping-like to field-like ratio GrjGi 
is influenced by the spin relaxation if it is anisotropic 
Xx 2 7 ^ Ayj. However, one sees that in the above equa¬ 
tion, the more detrimental factor is the penetration depth 
of the evanescent wave function = 1 / id+ + 9-) ^ 
h/y/2m{Vo — e) determined by the insulating gap Vb — e. 
The size of insulating gap in oxide insulators is typically 
of eV, yielding ^pen of the order of lattice constant, so 
it is unlikely that {Xx 2 TXy^} can be smaller than ^pen- 
In other words, unless the anisotropy of spin relaxation 
is associated with an energy scale (spin-orbit coupling, 
impurities, interface roughness, etc, which are beyond 
the scope of this approach) that is comparable to the 
insulating gap, the anisotropy is not crucial. With the 
assumption that ^pen Ayj} and a clean interface, 

from Eq. (4) one sees that the influence of {A 2 : 2 ,Ay 2 } on 
GrjGi is negligible. Physically, this means that the elec¬ 
tron can only penetrate EMI over a short distance within 
which its spin has yet started to relax. 


B. NM/FMM/oxide trilayer 


By matching the wave functions in Eq. (1), (II), and 
( 12 ) in the main text and their derivative in the two 
interfaces x = —b and a; = 0 , the scattering coefficients 
in NM/FMM/oxide trilayer are 
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where Uap, Zrrap-: and 7 ^ are those defined in Eq. (13) of 
the main text. 





























6 


II. DETAIL OF SPIN PUMPING 
CALCULATION 


other words, 


X2 = 



(a) 

5 II <<t)o II 



y2 = 




dt 


( 22 ) 


The infinitely small angle developed between S and 
(cr(a;,0))o after time lapse St is 




dS 

, 1 

dS 
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dt 
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(23) 


FIG. 3: (color online) (a) Schematics of the magnetization 
dynamics and the coordinate after time lapse St. (b) The 
magnetization dynamics induced by FMR, where the field¬ 
like component (green arrows) time-averages to zero, while 
the damping-like component (red arrows) does not. 


therefore 


lim sin0 = 

6t —^0 

lim cos 6 = 

it-s-O 


dS 

St = - 

dS 

dt 

1 . 
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dt 
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(24) 


The spin pumping can be calculated from our tunnel¬ 
ing formalism using adiabatic assumption, i.e., the mag¬ 
netization dynamics is much slower than any character¬ 
istic frequency scales in the problem. We start from the 
Bloch equation described by Eq. (9) in the main text and 
use NM/FMI bilayer as an example. When the magne¬ 
tization S is static, the conduction electrons leaked from 
NM into FMI are spin-split because of the exchange field 
rS • (T, causing an equilibrium spin density that can be 
described by Eq. (4) of the main text, where the incoming 
flux |Ap = Abi’l/.tol/a^ and the spin voltage at the inter¬ 
face /iQ are regarded as phenomenological parameters to 
ht the equilibrium spin density. The magnetization dy¬ 
namics modifies this equilibrium spin density as we now 
address. Suppose at f = 0, the conduction electron spin 
is in equilibrium with the magnetization, so 0 = 0 and 

(cr(x,0)) = (cr(x,0))o II S II Z 2 . (20) 

Because of the magnetization dynamics dS/dt, within an 
infinitely small time lapse St, the conduction electron spin 
is modified by 

St)) = 0))o + S{(t{x, St)) . (21) 


The modification of spin developed after the small time 
lapse St is, using Eq. (4) in the main text (ignoring spin 
relaxation as argued in Sec. 1.1) together with Eqs. (22) 
and (24), 
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So the time derivative in the Bloch equation is indepen¬ 
dent of St 
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If dS/dt is much slower than any characteristic frequency 
scale {e/h,Vo/h,T/h} in the problem, then one can as¬ 
sume that the system remains quasi-equilibrium so the 
wave functions we used are still applicable. This is gen¬ 
erally true in present experiments, since the magnetiza¬ 
tion dynamics, usually generated by ferromagnetic res¬ 
onance (FMR), is typically ~GHz, while the fermionic 
frequency scales {e/h,Vo/h, T/h} are usually above THz. 
Since dS/dt is pulling S away from {cr{x, 0))o, it is along 
X 2 , and S x dS/dt is along y 2 , as shown in Fig. 3 (a). In 


(cr^^)dS {<7^^)^ dS 
sin 9 dt sin 6 dt 


(26) 


Similarly, the other two terms in the Bloch equation are 


x{a) = ^ 

^ ‘^sd 


Tsf = 


Sjo') 

'^sf 


St 

'^Sf 


^ (gy^) 

dt sin 9 

dt sin 9 


dt sm 9 
dS (gy^) 


+ Sx 


dt sin 9 


(27) 
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which are proportional to St and hence vanish at dt —> 0. 
This also indicates that both the spin relaxation time Tgf 
and the s — d interaction time scale = h/TS are not 
detrimental to the spin pumping mechanism. The Bloch 
equation is then reduced to the continuity equation 

^ + dxjx = 0 . (at St —>■ 0) (28) 

By integrating over x and using j^o = 0, one obtains 
Eq. (10) of the main text. The spin pumping in the 
NM/FMM/oxide trilayer follows the same argument with 
the proper adaption of wave function and spin accumu¬ 
lation therein, and one arrives at the same Eq. (10). 
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